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1 Introduction 



One dimensional quantum spin systems have attracted much interest in recent years. On the 
theoretical side, soluble models and powerful methods such as bosonization combined with 
numerical results have been used to understand many aspects of the rich physics found in 
these systems. In particular the Bethe Ansatz solution of the S — 1/2 Heisenberg chain 
with nearest neighbour interaction has provided much insight into their properties. At the 
same time the fabrication of materials which are essentially quasi-one dimensional spin-1/2 
systems allows for experimental studies of these features. In addition to nearest neighbour 
exchange the effects of competing next-nearest neighbour interactions driving a spin Peierls 
transition ^-ff] and interchain interaction as in spin ladders f||-|8| have been considered to 
account for the properties of these compounds. Furthermore multi-spin exchange terms can be 
included without breaking the SU(2) symmetry of the system [pHTT[|. In general such terms are 
known to drive phase transitions in the Heisenberg chain opening a gap A between the singlet 
ground state and the lowest (triplet) excitation. 

New interesting properties can be observed when such systems are placed in a magnetic 



field ||i"2|, |T3|| : For H > A/gfig the system undergoes a transition from the gapped phase to an 
incommensurate phase with non-zero magnetization in which it has gapless excitations until 
the field is strong enough to lead to a saturated ferromagnetic state. Signatures of these 
transitions and properties of the gapless high-field phase have already been studied in several 
organic S = 1/2 spin Peierls and ladder systems |jl4|-|T6 . 



The subject of this paper is to study of the properties of this gapless phase. Our starting 
point is a soluble model of coupled spin-1/2 Heisenberg chains (0 < k < 1) ||T7H20[ 



N 

7~t — \ S2i-lS2i + S2i<SWlJ + J"i \S 2 iS 2 i+2 + S 2 i + iS 2 i + z 
i=l 

+J3 (S 2l x - S 2l (S 2i+1 x 5 2l+2 )) } (1.1) 



Ji = 2J(1 - k) , J 2 = Jk , J3 = 2J v / k(1 - k) . 

The three-spin exchange term guarantees integrability of the system for arbitrary ratio of 
antiferromagnetic pair exchange constants, i.e. < k < 1 (for more general cases see e.g. [f2~T| ). 
It breaks parity and time-reversal invariance, for a possible mechanism for the appearence of 



such a term see Ref . [3^] . Below we put J = 1 which fixes our scale of energy. In the following 
section we use the Bethe Ansatz solution to give a detailed account of the magnetic phase 
diagram and the critical properties of this model. Unlike the systems discussed above the 
model ( |1 . lp has gapless excitations for any k as long as the magnetization is not saturated. 



Furthermore, for sufficiently large k we find a transition into a second critical phase with 
different universality class at an intermediate magnetic field P0| |. Based on these findings we 
introduce a numerical method for the determination of the phase boundaries which allows to 
extend this discussion to nonintegrable generalizations of ( |1 . 1|) . Finally we apply this method 
to propose a phase diagram of the system with competing nearest and next-nearest neighbour 



exchange (i.e. without the three spin terms in fll.l)). 



2 The integrable model 

Starting from the ferromagnetic eigenstate with all spins up the spectrum of the integrable 
Hamiltonian ( |1 . 1| ) is obtained by adding magnons which are parametrized by the roots of the 
Bethe-Ansatz equations (BAE) JTj 



*3 ~ K + 2 \ / *3 + K + 2 \ TT Aj — Afc + 2 



with Ik = a/k/(1 — k). The corresponding state has magnetization A4 = N — M and energy 

M 2n 
E({\ j }) = 7iJ2( a i( X j + K) + ai(\ J -K)) , 7ra„.(A) = - 4A2 + ^ 2 . (2.2) 

i=i 

Generic solutions to the BAE are arranged in m-strings of complex rapidities A^ = x + i((m + 
l)/2 — j), j = 1, . . . , m with common real part. For the ground state only real solutions to (|2.1|) 
have to be taken into account. In the thermodynamic limit the ground state in a magnetic 
field H = h/gfiB is obtained by filling all states with negative dressed energy, defined in terms 
of the integral equation 



ei(A) + / dfi ci2(A — /i)e(yu) = h + n (ai(A + k) + ai(A — «)) . (2.3) 

Je{X)<0 

For h = one finds e\ < for all A, the corresponding ground state is parameterized uniquely 
by iV real rapidities, giving a singlet state with energy fi?(x) is the digamma function) 

/: -ln2- ( --ik) +^(l + ik) -* ( - + ik J 1 . (2.4) 



N 4 v ' V2 J v 7 V2 



The low-lying excitations are scattering states of spinons |23| and the low energy behaviour 
of the theory is that of a level-1 SU(2) Wess-Zumino-Witten model with central charge c = 1 
independent of k < 1 |2"4|. A small magnetic field breaks the SU(2) symmetry giving a c = 1 



Gaussian conformal field theory (CFT) with anomalous dimensions depending on h. Increasing 
the magnetic field h beyond 

, for < k < \ 

h>h c2 ={ I ~ " 4 2.5 

for i < k < 1 




one has ei(A) > resulting in a ferromagnetically polarized ground state [20|. 

For k > 1/4 an intermediate phase can be found for sufficiently strong fields h c \ < h < h c2 
which is characterized by 

ei(A) < for Ai < |A| < A 2 (2.6) 

(Ai 5 2 have to be determined as functions of the magnetic field from (|2.3|) .) Here the ground 
state is formed by filling two 'Fermi seas' of magnons with positive (negative) rapidities A. The 
presence of gapless excitations near the Fermi points at A = ±Ai, ±A 2 changes the low energy 
spectrum of the system which in this region has to be described in terms of two c = 1 Gaussian 
CFTs |20| very similar to a Luttinger liquid of electrons in one spatial dimension [2^, 26] . The 



asymptotic behaviour of correlation functions can be determined from the finite-size corrections 
to the energy of low lying states. Following Refs. |Z7],|28| we obtain (see Appendix [A]) 



2 

AE = ~QPP {Vl + lh) + JP + Ar) + V2{A * + A ^ (2 ' 7) 

where v± t2 are the Fermi velocities ([A.6Q of excitations near the points A\ )2 and 

± 1 f (x 12 AM 2 - x 22 AM 1 ) (z 12 AD 2 - ZggAgij 
1 ' 2{ 2(detX) T 2(detZ) 

± 1 f (x n AM 2 - x 21 AM x ) _ (z n AD 2 - z 21 AD 1 ) 



2 ' 2 \ 2(detX) + 2(detZ) J 1 ; 

are the conformal dimensions of the primary fields in the two Gaussian models expressed in 
terms of 2 x 2-matrices X and Z with elements ( |A.4|) and (|A.5|) respectively^ The integers 



AMi and ADi denote the difference of the quantum numbers ( |A.2[ ) between the excitation and 
the ground state. 

Taking the limit h \ h cl corresponding to A\ = the matrix elements of X and Z can 
be expressed in terms of the scalar dressed charge £ satisfying the integral equation £(A) = 
1 _ J A 1 2 rf /WA - V)t(v) as 

x 1 . U: J, ) • z : vl' ' • I < 29 > 




1 In other systems allowing for several types of massless excitations the relation X oc Z _1 allows to parame- 
terize the critical exponents in terms of a single dressed charge matrix |2£, 26,29,^0). Contrary to our claim in 
Rcf. pO] this is not possible in this model. 



This allows to simplify the expressions for the conformal dimensions ( |2.8| ) for h = h c \ + as 
Af = ]- {AMi ± [(1 - £(0))A£> 2 - ADi}} 2 



At = I C'-^^^ WI . (2.10) 



From ( |2.7| ) the correlation function of the operator with quantum numbers AMj, AD, is found 
to decay algebraically as with an exponent a = 2 (A^ + Aj[ + A^ + A^) p5| , [26|1 . Due 



to the nature of the ground state in the regime h c \ < h < h c2 several low energy processes 
contributing to C zz {x) = (S z (x)S z (0)) have to be considered (see Fig. |l|). The corresponding 
quantum numbers in (|2.8|) are 



(1) AL>i = AD 2 = 2 A Mi = AM 2 = 

(2) AL>i = -1 AD 2 = 1 A Mi = 1 AM 2 = 1 

1 ; (2.11) 

(3) ADi = 1 AD 2 = 1 A Mi = 1 AM 2 = 1 

(4) ADx = 2 AD 2 = A Mi = AM 2 = 

The leading asymptotic behaviour of C zz (x) is determined by the smallest exponent a. Ap- 
proaching the transition line h = /i c i(ac) from above we find that this exponent corresponds to 
the process (1) for 1/4 < /t < n c (h c i) and (3) for K c (h c i) < k < 1. At k = K c (h c i) ~ 0.61 or 
hdi^c) ~ 1.35 a crossover from an effective single chain to two chain behaviour takes place. 
Comparing the resulting exponents a to the one obtained in the single Gaussian theory for 
h / h cl we find from ( pA0| ) 



2£(A 2 ) 2 < h< h cl 

a=< 2(l-^(0)) 2 + 2e(A 2 ) 2 h\h cU K<n c (h cl ) (2.12) 

\ (1 + ^(A 2 )- 2 )(^(A 2 ) 2 + £(0) 2 ) h \ /i c i, « > K c (h cl ) 

This crossover in the long distance asymptotics of C zz [x) in the high field phase with h c \ < h < 
h c2 can be observed for any k > K c (h cl ): in the limit k — > 1 of decoupled chains the exponents 
due to the most relevant processes (1) and (3) in ( |2.11|) are = 1/£(A) 2 and = 2^(A) 2 



where A = (A 2 — Ai)/2. Again the dominant process at low energies is determined by the 
effective single-chain process with exponent c^ 3 -* for small fields with a crossover to the two- 
chain process (1) at a magnetic field to be determined from 2£(A) 4 = 1, or h ~ 1.29. Solving the 
integral equations numerically we obtain the complete phase diagram of the integrable chain 
in a magnetic field (Fig. 0). Depending on k and h the low energy properties of the system are 
those of a single or of two Gaussian CFTs. Increasing the magnetic field in the latter regime a 
crossover is observed corresponding a change from the intrachain process (3) to the interchain 



process (1) as the most important one at low energies. The typical dependence of the exponent 
a for fixed ratio of the exchange constants on the magnetic field is shown in Fig. |3|. 



3 Phase diagram of generalized spin ladders 



To investigate whether the magnetic phases found in the integrable system ( |1 . 1| ) persist when we 
relax the conditions on the pair and three spin exchange constants we consider the Hamiltonian 
(|1 . 1|) with general values of the exchange constants Jj. For magnetic fields strong enough to 
polarize the system completely (corresponding to h C 2 in the previous section) the excitations 
above the ferromagnetic ground state are spin waves. Their dispersion is easily obtained to be 



E± = J 2 cos k±yjjf cos 2 (A;/2) + ( J|/4) sin 2 k + const. (3.1) 

Below h C 2 the ferromagnetic state becomes unstable against creation of these magnons. De- 
pending on the coupling constants J, one finds either a single magnon mode at k = or two 
modes with wave numbers ±Q becoming soft at the transition into the paramagnetic phase — 
just as in the integrable model discussed in the previous section. The transition between these 
two scenarios occurs at the point where the low energy magnon modes (|3.1| ) have a quartic 
dispersion E- ~ k 4 for small wave numbers. This requires a ratio 

HH) 

for the pair exchange constants. For the integrable case (|1.1| ) this is the point (k, h) = (1/4, 3) 
in the phase diagram Fig. £2], i.e. the upper end point of the transition line between the two 
gapless phases. For a model without three spin exchange terms (J\, J2, -h) — (2(1 — k),k,0) 
Eq. ( |3.2|) implies a similar transition for k — 1/3 or J2/J1 = 1/4 (see also below). 

For magnetic fields h < h C 2 computation of spin wave dispersions as in ( [3.1|) is not suf- 
ficient due to strong quantum fluctuations in one spatial dimension. A possible criterion to 
determine the phase boundary at smaller fields is the field dependence of the zero temperature 
magnetization which shows a characteristic singularity on the transition line [pOfl . However, 
while this feature can be identified easily in the thermodynamic limit accessible for the Bethe 
Ansatz soluble model ( |1.1|) it cannot be used to analyze the numerical data due to the discrete 
set of magnetizations realized in a finite size system. Similarly, the dramatic changes in the 
low energy spectrum from a single Gaussian model to the form Q2.7| ) are difficult to see from 
the finite size data obtained by numerical diagonalization. 

As a possible method to identify the transition line h c i from numerical finite size data we 
propose the following method: having identified the phase for small ratios J2/J1 as similar in 



nature with the single Heisenberg chain J 2 = = J3 we compute the overlap 



0{J U J 2 , J 3 ) = \({Jx, 0, 0} I {J x , J 2 , J 3 }> I (3.3) 



between the finite size ground states | {Jx, J2, J3}) of ( |1.1| ) with fixed magnetization and the 
corresponding one of the Heisenberg chain. Clearly, this quantity will decrease from the value 
0(J\, 0, 0) = 1 for increasing next nearest neighbour and three spin interactions. In Fig. |] we 



present data for (|3.3| ) obtained for the integrable system fll.l ) with 16 spins in the sectors with 
different magnetization showing a sharp transition as a function of k. As shown in Fig. |2] these 
transitions provide an excellent numerical estimate of the critical field h c x- 

Away from the integrable line this sharp transition of the overlap (|3.3| ) as a function of 
J2I J\ is smoothed out. However, there are still well defined inflection points of the function 
0{Jx, J21 J 3) along lines in the space of parameters J 2j3 . This allows to estimate the position h c x 
of the transition in the nonintegrable model. In Fig. [5| the resulting phase diagram is shown for 
the system with nearest and next nearest neighbour pair exchange only, i.e. (Jx, J2, J3) = (2(1 — 
k), k, 0). Without an external magnetic field this model is known to have gapless excitations for 
< J2/J1 < 0.241 and in the limit of decoupled chains Jx — > [§-0]. For intermediate values 
of the next nearest neighbour exchange the system has a spin gap A leading to a plateau in at 
M z = in the magnetization curve extending to h = A. For strong magnetic fields h ~ h c2 the 
spin wave result ( |3.2|) shows a transition from a phase with one and two types of magnons with 
quadratic dispersion, respectively. For intermediate fields the K-dependence of 0(2(1 — k),k, 0) 
allows to locate the transition between these phases which is found to end near k ~ 0.5 at the 
transition into the gapped phase. 



4 Summary and Conclusion 

To summarize we have presented a detailed account of the magnetic phase diagram of an 
integrable model ( |1 . 1|) of generalized coupled spin chains. We have established two distinct zero 
temperature phases which can be described in terms of a single and two Gaussian conformal field 
theories respectively. The exact finite-size corrections in the spectra of low energy excitations 
in the latter phase allow to classify the possible critical exponents arising in the long distance 
asymptotics of two-point correlation functions. As an example the exponent in C zz (x) has 
been found to show a strong dependence on the magnetic field. In the same correlator we have 
found a crossover from an intrachain to an interchain process as the most relevant one at low 
energies. 



Finally we have extended our discussion to the possible magnetic phases of the system 
(|1.1|) with general exchange constants, in particular the spin chain with competing nearest 
and next-nearest pair exchange only. Combining known results on the nature of the zero field 
ground state with spin wave calculations and numerical simulations using the overlap ( [3 .3D as 
a criterion for the transition we have obtained the phase diagram of this system. Assuming 
that the phase found at large k is of similar nature as in the integrable model (this is certainly 
true in the large field limit) this gives strong constraints on a proper bosonization of the model 
in the gapless phase found for h > A. We should note, however, that by construction our 
analysis can not exclude the existence of additional phases at larger values of J%j J\ (which 
are not present in the integrable model). An study in the limit of two weakly coupled chains 
(Ji <C J 2 ) in a magnetic field may give additional insights on this region. 
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A Finite size spectrum of the integrable model 

To compute the finite-size spectrum of the integrable model in the regime h c i < h < h c2 we 
consider solutions of the BAE Q2.1J ) where all possible real solutions in the intervals Z_ = 
[A^~,Aj~] and X + = [A^,A^] with A^~ < A± are present (see Refs. [p7,|28|l). Their density 
satisfies the equation 



p(A) + jy + J jdfi a 2 (X - fi)p(ji) = a 1 (X + n) + <n(A - «) . (A.l) 

where j ± denotes integration over the interval T±. From p we define the quantum numbers 
(related to total magnetization and momentum) of the corresponding state as 

Mi r A * . , . . Di 



mi = — = J dpp(p) , di = — = [ I - I ) dpp(p) . (A.2) 



To express the finite size corrections to the energy of this state in terms of these quantum 
numbers, the matrix d{rrii, di}/d{A~j", A~} has to be computed. Defining function g(A|A) with 

g(X\A) + \ [ + f\dfi a 2 (X - p)g(p\A) = a 2 (X - A) (A.3) 



7 



and restricting ourselves to the case Af = ±A$ relevant for the Hamiltonian ( |1 . 1| ) we obtain 



1 dm,; 



X{j — 



and 



p(Ai) dAj 
1 ddj 



A, 



dfj, g{fj,\Aj) 



(A.4) 



Si4 - -IV 



-A, 



dfig(fj,\Aj) . 



Following Woynarovich (28[ these identities yield Q2.7|) and (|2.8| ) with the Fermi velocities 

1 dei (A) 



(A.5) 



-1 



(A.6) 



A=Aj 
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Figure 1: Low-energy processes contributing to C zz (x). 
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Figure 2: Magnetic phase diagram of the integrable model. The dashed line marks the crossover 
between the intrachain and interchain process as dominant low energy contribution to C zz , the 
symbols denote the position of the phase transition as determined from numerical diagonaliza- 
tion of systems with 16 and 20 spins using the singularities in the operlaps ( |3.3| ) . 
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Figure 3: Critical exponent a of C zz [x) for k = 3/4 as a function of the magnetic field. The 
cusp at h ~ 1.31 is a consequence of the crossover from the intrachain (3) to the interchain 
process (1) in Fig. |I[ 




Figure 4: Overlap matrix elements ( |3.3|) for the integrable model with 16 spins as a function of 
k: the data at small k correspond to states with magnetization M z = 6, 4, 2, (top to bottom), 
the lines are guides to the eye. 
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Figure 5: Phase diagram of the spin chain model with nearest and next-nearest neighbour pair 
exchange 2J 2 / J\ = k/(1 — k): the upper line denotes the saturation field h c2 obtained from spin 
wave calculation, the lower one denotes the transition at h — A into the dimerized phase with 
a spin gap (0.325 < k < 1 for h — > 0). Squares (diamonds) denote the transition determined 
numerically from the overlap criterion ( |3.3| ) for systems with 20 (16) spins. 
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